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Abstract 

The order a mixed QCD and weak corrections to top quark pair production by quark antiquark anni- 
hilation are computed, keeping the full dependence on the t and i spins. We determine the contributions 
to the cross section and to single and double top spin asymmetries at the parton level. These results are 
necessary ingredients for precise standard model predictions of top quark observables, in particular of 
top spin-induced parity-violating angular correlations and asymmetries at hadron colliders. 
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One promising tool for investigating the so-far relatively unexplored dynamics of top quark production 
and decay, once high statistics samples of t and/or t quarks are available, are observables associated 
with the spins of these quarks. As far as QCD-induced tt production and decay at hadron colliders is 
concerned, theoretical predictions for differential distributions including the full dependence on the f,f 
spins are available at NLO in the QCD coupling [1,2]. 

For full exploration of sizeable, respectively large tt data samples that are expected at the Tevatron and 
at the LHC the standard model (SM) predictions should be as precise as possible. Specifically weak 
interaction contributions to tt production should be taken into account. Although they are nominally 
subdominant with respect to the QCD contributions they can become important at large tt invariant mass 
due to large Sudakov logarithms (for reviews and references, see e.g. [3,4]). 

SM weak interaction effects in hadronic production of heavy quark pairs were considered previously. The 
parity-even and parity-odd order a^a vertex corrections^ were determined in [5] and in [7], respectively 
(see also [6]). In ref. [7] also parity-violating non-SM effects were analysed. The box contributions to 
qq tt were not taken into account in these papers. In [8] the weak contributions to the hadronic bb 
cross section, including these box contributions, were computed. 

In this Letter we report on the calculation of the mixed QCD and weak radiative corrections of order 
a^a to the (differential) cross section of tt production by quark-antiquark annihilation, keeping the full 
information on the spin state of the tt system. These results are necessary ingredients for definite SM 
predictions, in particular of parity-violating observables associated with the spin of the (anti)top quark. 
In the following we first give some details of our calculation. Then we present numerical results for the 
cross section and for several single spin and spin-spin correlation observables. 

Top quark pair production both at the Tevatron and at the LHC is dominated by the QCD contributions 
to qq tt and gg tt, which are known to order a^. Due to color conservation there are no a^a Bom 
level contributions to these processes. The leading contributions involving electroweak interactions are 
the order Bom terms for qq tt and the mixed contributions of order a^a. For the quark-antiquark 
annihilation processes, which we analyze in the following, this amounts to studying the reactions 

(l{P\)+^{P2) -^t{ki,St)+t{k2,Si), (1) 

l{Pi) +q{P2) t{ki,st) +t{k2,St) +gik3), (2) 

Here pi, p2, k\, ki, and ^3 denote the parton momenta. The vectors St, sj, with sj = = —\ and 
ki- St = k2- sj = describe the spin of the top and antitop quarks. All quarks but the top quark are taken 
to be massless. 

The respective contributions to the differential cross section of O are of the form 

\M2{p, k, St, St) p + a^a hM2{p, k, st^sj) , (3) 

where M2 corresponds to the y and Z exchange diagrams. As we are interested in this Letter in particular 
in parity-violating effects, we take into account only the mixed QCD and weak contributions to dM2 
and to © in the following. The photonic contributions form a gauge invariant set and can be straight- 
forwardly obtained separately. The contributions to dM2 are the order two-gluon box diagrams in- 
terfering with the Bom Z-exchange diagram, and the Z gluon (g) box diagrams and the diagrams with 

'Here and a denote the strong and electromagnetic couplings, and the weak coupling is = a/sin^Qw 
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the weak corrections to the qqg and gtt vertices interfering with the Born gluon exchange diagram. The 
ultraviolet divergences in the vertex corrections are removed using the on-shell scheme for defining the 
wave function renormalizations of the quarks and the top quark mass nit. 

The respective contributions to the differential cross section of © are of the form a?^ahM-i{p,k,St,St) 
and result from the interference of the order with the order gse^ gluon bremsstrahlung diagrams. 
The box diagram contributions to © contain infrared divergences due to virtual soft gluons. They are 
canceled against terms from soft gluon bremsstrahlung. As as consequence of color conservation both 
the sum of the box diagram contributions to and hM^, are free of coUinear divergences. 
We have extracted the IR divergences, using dimensional regularization, with two different methods: a 
phase space slicing procedure (as in [2]) and, alternatively as a check, we have constructed subtraction 
terms that render the three particle phase space integral over the subtracted term [^M^]subtr finite. When 
calculating observables, in particular those given below, both methods led to results which numerically 
agree to high precision. 

We have determined © and 6fM'3, respectively their infrared-finite counterparts, analytically for arbri- 
trary t and t spin states. From these expressions one may extract the respective production spin density 
matrices. These matrices, combined with the decay density matrices describing semi- and non-leptonic t 
and t decay [9] then yield, in the tt leading pole or narrow width approximation, standard model predic- 
tions for distributions of the reactions qq^tt^ bb + 4f {+g) (/ = q.i.Vg) with the t and f spin degrees 
of freedom fully taken into account. 

The expressions for © and dMj, are rather lengthy when the full dependence on the t and t spins is kept, 
and we do not reproduce them here. We represent these contributions to the partonic cross sections and 
to several single and double spin asymmetries, which we believe are of interest to phenomenology, in 
terms of dimensionless scaling functions depending on the kinematic variable r\ = — 1, where s is 

the qq cm. energy squared. The inclusive, spin-summed qq cross sections for ([l]), © may be written, to 
NLO in the SM couplings, in the form 

o^-oS^^VSaf + 8<,, (4) 

where the first and second term are the LO (order a^) and NLO (order a,) QCD contributions [10, 1 1], 
[12], and the third term is generated by the electroweak contributions © and SfWs described above. We 
decompose this term as follows: 

5a,^-(f,m2) = ^[a/^^^Ti) -fa2/ii^(Tl)] • (5) 

We have numerically evaluated the scaling functions (rj) - and those given below - and parameterized 
them in terms of fits which allow for a quick use in applications. In the following we use = 178 GeV, 
mz = 91.188 GeV, and sirP'Qyv =0.231. In Figure 1 and Figs. 3-9 below we use nin = 1 14 GeV for the 
mass of the standard model Higgs boson. The dependence on the Higgs boson mass is shown in Fig. 2 

in the case of /j^ for two values of mn [13]. 

In Fig. [T]the functions f^l^ are displayed as functions of T] for annihilation of initial massless partons 
qq of the first and second generation with weak isospin ±1/2. As expected the a^a corrections are 
significantly larger than the the lowest order photon and Z boson exchange contributions. The correction 
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dll) to the qq cross section has recently been computed also by [14]. We have compared our results and 
find excellent numerical agreement. 

In Fig. |2lthe contributions to /j^"* of the initial and final vertex corrections and of the box plus gluon radi- 
ation terms are shown. These two contributions are separately infrared- finite. This figure clearly shows 
that the latter contributions should not be neglected. This statement holds also for the spin observables 
discussed below. We have numerically compared the contributions of the initial and final vertex correc- 
tions to ^ relative to the order QCD Bom cross section with the results Figs. 9 and 10 of [5] and find 
agreement. Fig. |21also shows that for T] < 10 the dependence on the Higgs boson mass is significant. 
The contribution of the corrections shown in Fig. [l]to the tt cross section at the Tevatron is very small. 
This is mainly due to the fact that the order a^a corrections change sign for larger Tj. The significance 
of the contributions can be enhanced by suitable cuts, e.g. in the tt invariant mass. 
In Fig. |3the order a^, a^, and the a^a contributions to the cross section (|H) are shown as functions of 
r|. In these plots as{mt) = 0.095 and a(mz) = 0.008 was chosen. One sees that for t] > 1 the mixed 
corrections become of the same size or larger in magnitude than the NLO QCD contributions, and at 
T] ~ 10 the a^a contributions are already about 15 percent of the LO QCD cross section. These regions 
can be investigated by studying the distribution of the tt invariant mass Mtt- A value of, say, T] ~ 10 
corresponds roughly to Mjt ~ 1 TeV. For a quantitative discussion the reaction gg tiX must, of course, 
also be taken into account. At the LHC, where such studies may be feasible, this is the dominant channel. 

Next we consider spin asymmetries. Denoting the top spin operator by St and its projection onto an 
arbitrary unit axis a by S/ • a we can express its unnormalized partonic expectation value, which we 
denote by double brackets, in terms of the difference between the "spin up" and "spin down" cross 
sections: 

2((Sra)), = o,-(T)-a,-(i). (6) 

Here i = qq and the arrows refer to the spin state of the top quark with respect to a. An analogous formula 
holds for the antitop quark. It is these expressions that enter the predictions for (anti)proton collisions. 
There are two types of single spin asymmetries ® : parity-even, T-odd asymmetries and parity-violating, 
T-even ones. The asymmetry associated with the projection S/ onto the normal of the q, t scattering plane 
belongs to the first class. It is induced by the absorptive part of d9\{2, but also by the absorptive part of 
the NLO QCD amplitude. (This was calculated in [15, 16].) The weak contribution is even smaller than 
the one from QCD which is of the order of a few percent. For the sake of brevity we do not display it 
here. 

The P-odd, T-even single spin asymmetries correspond to projections of the top spin onto a polar vector, 
in particular onto an axis that lies in the scattering plane. Needless to say, they cannot be generated 
within QCD; the leading SM contributions to these asymmetries are the parity-violating pieces of eq. 
^ and 5fAf3 above. We consider top spin projections onto the beam axis (which is relevant for the 
Tevatron), onto the helicity axis (of relevance for the LHC), and for completeness also onto the so-called 
off-diagonal axis, which was constructed to maximize tt spin correlations in the the qq channel [17]. 
Naively, one might define these axes in the cm. frame of the initial partons. However, the observables 
St ■ a are then not coUinear-safe. (The problem shows up once second-order QCD corrections are taken 
into account.) A convenient frame with respect to which coUinear-safe spin projections can be defined 



3 



is the tt zero-momentum-frame (ZMF) [2]. With respect to this frame we define the axes 



a = b = p, (beam basis), (7) 
a = b = d, (off-diagonal basis), (8) 
a = — b = k, (helicity axes), (9) 

where k denotes the direction of flight of the top quark in the tt-ZMF and p is the direction of flight of 
one of the colliding hadrons in that frame. The direction of the hadron beam can be identified to a very 
good approximation with the direction of flight of one of the initial partons. The unit vectors b are used 
for the projections of the spin of the t quark. The vector d is given by 

, -p+(l-T)(p.k)k^ (10) 



^ l-(p-k)2(l-y2) 

where E /m. 

The unnormalized expectation values of St • a are again conveniently expressed by scaling functions 

((2S. -a))^, = ^[a4°/)(Ti) +a?4^'")(il)] . (11) 

The results for the scaling functions corresponding to the three axes above are shown in Figs. For 
the beam and off-diagonal axes the asymmetries are almost equal both at LO and at NLO, up to sign, 
and for weak isospin % = —1/2 quarks the a^a corrections are significantly larger than the LO values. 
This is in contrast to the helicity basis where the LO and NLO terms shown in Fig. |6lare of the same 
order of magnitude. Moreover, in this basis the LO and NLO terms cancel each other to a large extent 
for % = — 1 /2 quarks in the initial state. Because tt production at the Tevatron occurs predominantly by 
qq annihilation, these results imply that, when it comes to searching for SM-induced parity-violating t 
or t spin effects, the reference axes of choice should be p or d. The SM effect is quite small; for suitably 
chosen Mfj mass bins one gets asymmetries of the order of 2 percent. This leaves a large margin in the 
search for new physics contributions. 

Finally we consider top-antitop spin-spin correlations. For the sake of brevity we concentrate here on 
parity- and T-even ones which are generated already to lowest order QCD. For the Tevatron these spin 
correlations (including NLO corrections) are largest with respect to the beam and off-diagonal bases, 
while for the LHC the helicity basis is a good choice^. In addition, a good measure for the tt spin 
correlations at the LHC is, in the case of the dilepton channels, the distribution of the opening angle 
between the charged leptons from semileptonic t and t decay. At the level of tt this amounts to the 
correlation S/ ■ Sf (for details, see [2]). Therefore we consider the following set of observables: 

Oi=4(p.Sf)(p-S,-), (12) 

02 = 4(d-S,)(d-S,-), (13) 
03 = -4(k-S,)(k-S,-), (14) 



For the LHC, a basis has been constructed [18] which gives a QCD effect which is somewhat larger than using the helicity 
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04 = 4S,-Sr, (15) 

where the axes are as defined in eqs. Q - (|9l) in the tt ZMF and the factor 4 is conventional. The unnor- 
malized expectation values of these observables correspond to unnormalized double spin asymmetries, 
i.e., to the following combination of spin-dependent cross sections: 

((0/,)), = o,(TT)+a,(ii)-o,(Ti)-o,-(iT), (16) 

and here i = qq. The arrows on the right-hand side refer to the spin state of the top and antitop quarks 
with respect to the reference axes a and b. 

Again we compute the and weak-QCD contributions of order a^a to (fT6b and express them in terms 
of scaling functions: 

These functions are plotted in Figures |7] - |9l respectively. As is the case in QCD, in qq annihilation 
the spin correlations in the beam and off-diagonal bases, and in the projection (fTSt are not very much 
different from each other. The size of the mixed corrections dTTl is typically only a few percent compared 
with the QCD contributions [2] to (fT6t . 

The unnormalized expectation value of O4 given in eq. dTSt is equal to the respective contribution 5o^ 
displayed in Fig. [H The reason is as follows. The S matrix elements for the reactions qq tt{g), to 
the order in the couplings considered above, contain all possible partial wave amplitudes. However, in 
the expectation value of the parity-even operator dTSt only the ^Si (in the Born terms and the terms of 
order a^a) and ^Pi components (in the Born term from Z boson exchange) of the tt state contribute, as 
a closer inspection shows. It is then a simple exercise to show that the normalized expectation value 
of O4 is equal to one in this case. Thus, its unnormalized expectation value is equal to 5o^, which is 
confirmed by explicit calculation. (See [19] for similar considerations.) 

Another interesting class of asymmetries are parity-violating double spin asymmetries of the form 
5A(a,b) = 0/(ti) — 0/(it)- They are generated by the parity- violating pieces of © and of hM^, above. 
In [7] an observable of this form was computed in the t,t helicity basis within the SM, with box plus 
gluon contributions not taken into account, and in some SM extensions. 

In addition, the absorptive parts of lead to T-odd spin-spin correlations, both P-even and odd ones. 
These are, however, very small effects, and we do not display them here. 

The single and double spin asymmetries are reflected in respective angular distributions/asymmetries of 
the t and i decay products. In particular they contribute to the one- and two-particle inclusive decay 
distributions o~^Ja/Jcos0i and ^ J a/ (J cos 0i J cos 62), where 9i,02 are the angles between the 
direction of flight of a ? and t decay product, respectively, and a chosen reference direction. Suitable 
reference directions are the axes introduced above. The effects are largest if the charged lepton(s) from 
t and/or t decay is (are) used as spin analyzer(s). The weak contributions can be enhanced with respect 
to the pure QCD effects by suitable cuts in the tt invariant mass. Numerical studies, including the weak 
contributions to gg — »• tt will be given elsewhere [20] . 

In summary we have computed the mixed QCD and weak corrections to top quark pair production by 
quark antiquark annihilation, keeping the full dependence on the t and t spins. These results, combined 
with our previous QCD results and with the mixed contributions to gg tt, will allow for detailed pre- 
dictions of top quark observables, in particular of top spin-induced angular correlations and asymmetries 



5 



within the standard model [20]. Specifically, the results of this letter are necessary ingredients for SM 
predictions of parity-violating observables associated with the spin of the (anti)top quark. We believe 
that such observables, which are very sensitive to non-SM parity-violating top quark interactions, will 
become important analysis tools once sufficiently large tt data samples will have been collected. 
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Figure 1: Dimensionless scaling functions f^!l\v[) (dashed), /IV (t]) (solid) that determine the parton 
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qq 



cross section ^ for q = d type quarks. The dash-dotted and dotted lines correspond to the respective 
functions for = m type quarks. The Higgs boson mass is put to 1 14 GeV. 
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Figure 2: Contributions of the initial and final vertex corrections for m//= 1 14 GeV (dashed), for mn = 
250 GeV (dash-dotted), and of the box plus gluon radiation terms (solid) to /jj (rj) . 
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Figure 3: Contributions of the LO (solid) and NLO QCD (dashed) contributions (taken from [12]) and 
of the mixed a^a contributions (dotted and dash-dotted line refers to initial d-type and u-type quarks, 
respectively) to the cross section Q in units of 1/m^, and niH = 1 14 GeV. 
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Figure 4: Scaling functions h!"^q^\r\) (dashed), h)^q^\r\) (solid) that determine the expectation value (fTTT) 
for the beam axis in the case ofq = d type quarks. The dash-dotted and dotted lines correspond to the 
respective functions for ^ = m type quarks, nin = 114 GeV. 
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Figure 5: Scaling functions hg^'''{r\) (dashed), h^^^" (solid) that determine the expectation value (fTTT) 
for the off-diagonal axis in the case of q = d type quarks. The dash-dotted and dotted lines correspond 
to the respective functions for q = u type quarks, m// = 1 14 GeV. 
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Figure 6: Scaling functions hg^'''{r\) (dashed), Z?^^'" (t]) (solid) that determine the expectation value (fTTl) 
for the helicity axis in the case of q = d type quarks. The dash-dotted and dotted lines correspond to the 
respective functions for ^ = m type quarks, m^j = 1 14GeV. 
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Figure 7: Scaling functions g^^-' ' (rj) (dashed), (rj) (solid) that determine the expectation value (fTTl) 
for the beam basis in the case of q = d type quarks. The dash-dotted and dotted lines correspond to the 
respective functions for ^ = m type quarks, mn = 114 GeV. 
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Figure 8: Scaling functions g"^^"^ (rj) (dashed), (solid) that determine the expectation value (fTTl) 



for the off-diagonal basis in the case of q = d type quarks. The dash-dotted and dotted lines correspond 
to the respective functions for = m type quarks, nin = 114 GeV. 
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Figure 9: Scaling functions g[^q^\v[) (dashed), g^qq^\v[) (solid) that determine the expectation value (fTTl) 
for the helicity basis in the case of q = d type quarks. The dash-dotted and dotted lines correspond to 
the respective functions for ^ = a type quarks, mh = 114 GeV. 
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